Quantum law of rare events for systems with Bose-Einstein statistics 



D. Sokolovski 1 ' 2 

1 Departmento de Quimica-Fisica, Universidad del Pais Vasco, UPV/EHU, Leioa, Spain 
2 IKERBASQUE, Basque Foundation for Science, E- 480 11 Bilbao, Spain 
(Dated: December 18, 2012) 

In classical physics the joint probability of a number of individually rare independent events is 
given by the Poisson distribution. It describes, for example, unidirectional transfer of population 
between the densely and sparsely populated states of a classical two-state system. We derive a 
quantum version of the law for a large number of non-interacting systems (particles) obeying Bose- 
Einstein statistics. The classical low is significantly modified by quantum interference, which allows, 
among other effects, for the counter flow of particles back into the densely populated state. Suggested 
observation of this classically forbidden counter flow effect can be achieved with modern laser-based 
techniques used for manipulating and trapping of cold atoms. 
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In classical physics and statistics, probability for a 
number of individually rare events is universally given 
by the Poisson distribution (see, for instance, [1 j). For 
example, it is obeyed by a classical gas escaping into an 
empty space through a penetrable membrane. With the 
number of atoms N large, and the transition probabil- 
ity made proportionally small, the number of escaped 
atoms is governed by the Poisson law, with the number 
of atoms recaptured into the original reservoir vanishing 
as N — >• oo. The validity of the Poisson distribution de- 
pends on that one can, in principle, know not only how 
many but also which of the atoms have escaped. Quan- 
tum mechanics offers a different possibility: for identical 
particles one is allowed to know only the number of the 
escapees, and not their identities. While it is well known 
that both Fermi-Dirac and Bose-Einstein symmetries of 
a wave function may lead to non-poissonian effects in the 
full counting statistics of otherwise independent particles 
[2] -[8], the failure of the Poisson law in the limit of rare 
events is less obvious. The subject of this Letter is the 
general question of what replaces the classical Poisson 
law in a quantum situation where only the total number 
of rare events, but not their individual details, can be 
observed?. 

We specify to the case of many non-interacting bosons, 
each of which may occupy one of the two available states. 
Such systems are also of practical interest, e.g., for their 
potential applications as detectors. For example, if the 
transmission amplitude between two connected cavities 
is influenced by a passing particle, the change observed 
in the photonic current would announce the particle's 
arrival. In a similar way, atomic current of a weakly in- 
teracting Bose-Einstein condensate (BEC) trapped in a 
double- or multi-well potential (see Fig.l) can be used to 
gain information about the state of a qubit coupled to 
the BEC [9] -[11]. Detailed analysis of the work of such 
hybrid bosonic devices must take into account in which 
manner, and how frequently, the bosonic sub- system is 
observed, and will be given elsewhere. 
We note that the problem is fundamentally different from 
that of the frequently studied coined quantum walk [12] , 
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FIG. 1. (Color online) Double- well trap containing N atoms. 
The central barrier is lowered to allow tunnelling between the 
states |1) and |2). 

where interference between virtual paths available to a 
single particle modifies the classical Gaussian distribu- 
tion. In our case, modification of the classical law is a 
many-body effect, specific to the Bose-Einstein statistics. 

We start by constructing transition amplitudes for a 
single particle, which can occupy one of the two levels 
in an asymmetric double well potential. In terms of the 
Pauli's matrices, the Hamiltonian reads 

H = ea z + £a x + r]a y , (1) 

where the spin states |1) and |2), aligned up and down 
the z-axis, represent an atom in the right and left state, 
respectively, e is the difference between the energies of 
the states, and £ and r] together define the tunnelling 
matrix element, T = £ + irj. For the evolution operator 
U(t) = exp(-iHt) we have 

U(t) = Icos{ujt) — iuj~ 1 [ecr z + £a x + rja y ] sin(cjt), (2) 

with oj = y^e 2 + £ 2 + t? 2 , and its matrix elements are con- 
veniently written as (a star indicates complex conjugate) 

Un = y/l -pexp(ia) = Z7 2 * 2 (3) 
U12 = ^/pexp(if3) = -UZn 
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FIG. 2. A diagram showing the region of summation 
in Eqs. Q and (|9| (filled dots). Each dot contributes 
C™C?- m p u+fl (l -py*-^-" for distinguishable particles, and 



CG^^t-l)^^^""^"""" f o r identical bosons. 



where, with our choice of the basis, p(t) = (£ 2 + 
rj 2 ) sm 2 (uot) / uj 2 is one-particle transition probability, 
a(t) = — arctan[e tan(cj£)/u;] and f3 = —tt/2. 

For a total of N particles, we wish to evaluate the tran- 
sition probabilities Pm f< _ m (i) for starting with m parti- 
cles in the state |1) and ending, after a time £, with m! 
particles in the same state. It is instructive to begin with 
a brief discussion of the case where all particles are con- 
sidered distinguishable. The problem is equivalent to a 
classical TV-coin one: given that each coin changes its 
state with a probability p, and m coins initially oriented 
heads up, what is the probability to have m! heads up 
after each coin has been tossed once? The result can 
be achieved by moving v coins from tails to heads, and \i 
coins from heads to tails, provided v — ji = m! — m. Sum- 
ming the corresponding probabilities, while taking into 
account the number of ways to choose the coins which 
change their state, yields 
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is the binomial coefficient, and <L 



J l = l\(k-l) 

is the Kronecker delta. Depending on 7V,ra and ra', the 
sum m Eq.Q may contain a different number of terms, 
corresponding to the number of 'pathways' connecting 
the initial and final states (filled dots in Fig. 2). In the 
rare events (RE) limit 
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(5) 



it is sufficient to retain only the leading \i = terms (the 
lowest dot in the diagram in Fig. 2) in Q. Using the rela- 
tion UrriN^ooCm = N m /m\ for the remaining binomial 
coefficient, yields the expected Poisson distribution, 
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FIG. 3. (Color online) Poissonian probability p^'^m f° r dis- 
tinguishable particles, as given by Eq.Q, for N = 10 5 and 
w = 3. The white line marks m — m . 



with intuitively appealing properties. Indeed, in the case 
of a symmetric trap, e = 0, reducing the transition proba- 
bility in ([5| also makes the Rabi period 2tt/uj after which 
the system must return to its initial state, extremely 
large. Now, for allt << 2ty/uj : the evolution can be con- 
sidered approximately irreversible, with number of parti- 
cles q escaping into the right trap independent of number 
of particles, m, already there. Low probability of each in- 
dividual event, and much lower population in the right 
well make re-crossings from right to left extemely unlikely 
(see Fig. 3). In particular, after detecting m particles in 
the right well, one never finds it empty again, as the 
probability Po^ m (the left upper corner in the diagram 
in Fig. 2), vanishes as (w/N) m exp(— w). One might ex- 
pect a similar argument to be also valid should distinct 
particles be replaced with non- interacting bosons. Next 
we will show that this is not the case. 

For identical bosons we have a quantum version of the 
iV-coin problem: after a toss each coin changes its state 
from \i) to z, j = 1,2 with the probability amplitude 
Uji, and we must sum amplitudes rather than probabili- 
ties over all pathways leading to the same final state. The 
state of the system with any m coins displaying heads is 
given by a symmetrised wave function 
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(7) 



where i = 1, 2 denotes the state of the j-th particle, 
and the sum is over different ways to ascribe to m 
of the N indices ij the value of 2, and to the remaining 
N — m ones the value of 1. After all coins are tossed once 
each individual term in the sum of Eq.([7|) contributes to 
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with the region of summation illustrated in Fig. 2. Since 
Eq.(7|) contains such terms, the probability to have 
m! heads up after the toss is P^ m = (C ? ^) 2 |/(m / <— 
m, 7V)| 2 , m,m' = 0, 1,...N, which, with the help of 
Eq.fpl), can be expressed in terms of the Jacobi poly- 
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(9) 



As in the case of distinguishable particles, P^^ m de- 
pends only on the one-particle transition probability p, 
and not of the phases a and j3 of the matrix elements of 
of Uij in Eq.|3| [16]. We note also that in the special case 
of tunnelling into an initially empty well, m = 0, there is 
only one pathway (moving exactly m particles from left 
to right), and transition probabilities for distinguishable 
particles and identical bosons coincide, 
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as was pointed out earlier in the Refs. [10] and [TT] . 
More interesting, however, are the transitions affected 
by the interference effects which, as we will demonstrate, 
persist even in the RE limit Indeed, since the sum in 
Eq.([8| contains ^/p rather than p, the restriction to only 
\± = (T terms is no longer justified. Thus, after taking the 
limit ([5]), we have (q = m' — m) 
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Equation (11), which is our central result [17 , replaces 
the classical Poisson law ([6| for non-interacting identical 
bosons. Some of its properties are counterintuitive, as 
is shown in Fig. 4. Firstly, unlike the Poisson distribu- 
tion (pi, Pm>m °f Eq.(ll) is highly structured, as a result 



of the interference between the pathways. Secondly, it 
allows for total or partial recapture of the few particles 
initially held in the right well back into the densely pop- 
ulated left well, contrary to the simple argument based 
on improbability of such an event. 

The probability for all m bosons to cross into the left 



well, P(^L m , contains only one term in the sum ( 11 ) [(/i 



0) in the diagram in Fig. 1], 

P P l m = w m ex V (-w)/m\. (12) 
As a function of the number of recaptured atoms m, 
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FIG. 4. (Color online) Non-poissonian probability P r ^/^ m for 
identical bosons, as given by Eq.|9]), for N — 10 5 and w = 3. 
The section of the surface indicate by the arrow corresponds 
also shown in Fig. 4. 
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FIG. 5. (Color online) The probability P ^ m in Eq.^ that 
all m particles initially in the right well, m « N, N » 1, 
would cross to the left, leaving the right well empty. The solid 
lines are the corresponding Poisson distributions (12) 



it is a Poisson distribution shown in Fig. 5 (apparently 
so ubiquitous that after having been evicted from one 
part of this paper, it immediately reappears in another, 
albeit in a different context). The recapture process ex- 
hibits certain resonance-like behaviour. The number of 
particles most likely to be readmitted to the left well, 
m « w, equals the mean number of distinguishable par- 
ticles crossing into the right well under the same con- 
ditions. For m > 0, there are two interfering scenarios 
leading to only one particle being left in the right well 
[points (/i = m — 1, v = 0) and (/i = m, v = 1) in Fig. 2], 
and the corresponding probability is bimodal as shown 
in Fig. 4. Similarly, the probability P^^ m to retain the 
same number of atoms in the right well builds up from 
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FIG. 6. (Color online) The probabilities Pi^ m to end up with 
just one atom in the right well (circles, long dashed), and 

(triangles, dashed), as functions of m, 



Using the available laser technology [3] one can create 
quasi-one-dimensional box-trap with two strong endcap 
lasers providing the potential walls shown in Fig.l. The 
box is divided in two by adding a third laser in the mid- 
dle, and the left well is populated with a large number TV 
of weakly interacting atoms, while, say, three atoms are 
introduced into the right well. Following this, the middle 
laser beam is slightly weakened to allow transfer of atoms 
between the wells. It is restored after a time r, such that 
p(r) ~ £ 2 r 2 « 3/7V, and the number of atoms in the 
right well is measured, e.g., by a technique described in 
[T8] . Then, no matter how large is TV, the well is found 
empty with the probability P(^L m ~ 0.23 (c.f., Fig. 5), 
i.e., in just under a quarter of all cases. 

In summary, it is shown that for non- or weakly in- 
teracting bosons quantum interference between different 
scenarios leading to the same final state modifies the clas- 
sical Poisson law of rare events, and leads to significant 
observable effects not present in classical statistics. 



m + 1 interfering terms and also shows an oscillatory pat- 
tern (see Fig. 6). 

To conclude, we suggest a simple experimental setup to 
test the re-capture property of the bosonic distribution. 
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